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^ ■ Abstract 

o 

f-| ' The hyper fine sphttings in heavy quarkonia are studied in a model- 

QhI independent way using the experimental data on di-electron widths. 

o ' Relativistic correlations are taken into account together with the smear- 

1^ ■ ing of the spin-spin interaction. The radius of smearing is fixed by the 

known J/ip — r]c{lS) and tp{2S) — r]'^{2S) splittings and appears to be 

small, rss — 0.06 fm. Nevertheless, even with such a small radius an 

^ ■ essential suppression of the hyperfine splittings (~ 50%) is observed 



in bottomonium. For the nS bb states (n = 1, 2, . . . , 6) we predict the 
values (in MeV) 28, 12, 10, 6, 6, and 3, respectively. For the 3S and 
45 charmonium states the splittings 16(2) MeV and 12(4) MeV are 
obtained. 



1 Introduction 

At present two spin-singlet S-wave states ric{lS) and r]c{2S) have been dis- 
covered P-[SI- Still, no spin-singlet r]b{nS) levels have been seen jlj. The- 
oretically the masses of the rj^^nS) were predicted in many papers |3]-[TT], 

*E-mail:badalian@itep.ru 



however, the calculated hyperfine (HF) splittings, 

AnpinS) = M{n^Si) - M{n^So), (1) 

vary in a wide range: from 35 MeV up to 100 MeV for the bb IS state and 
for the 2S state between 19 MeV and 44 MeV [TT] . 

However, modern theoretical treatments and current experiments taken 
together, have produced well-established limits on the factors which deter- 
mine the spin-spin potential Vnplr) in heavy quarkonia. First of all, the wave 
function (w.f.) at the origin for a given n^Si{cc or bb) state can be extracted 
from the di-electron width which, as well as some ratios of leptonic widths, 
are now measured with high accuracy ^21 1^ • The quark masses, the pole 
(current) mass, present in the correct relativistic approach, and the con- 
stituent mass, used in the nonrelativistic or in more refined approximations, 
is also known with good accuracy P^ ll5j. Therefore the only uncertainties 
comes from two sources. 

First, in perturbative QCD a strict prescription as to how to choose the 
renormalization scale /i in the strong coupling anp, which enters Vhf(^), is 
not well established for a bound state, especially for higher excitations. 

Secondly, a smearing of the spin-spin interaction is considered to be due 
to relativistic effects [TU] but the true size of the smearing radius r^g is still 
not fixed. 

In our calculations the radius r^s is taken to fit the J/ip — Vci^^) s^nd 
i/j{2S) —rjc{2S) splittings. We show that to reach agreement with experiment 
the smearing radius should be Tgs < 0.06 fm. Our value Tgs = 0.057 fm 
practically coincides with the number used in Ref. ^U]. However, in spite of 
this coincidence the splitting Ai = T(IS') — rih{lS) = 28 MeV found in our 
calculations appears to be two times smaller than that in Ref. |iOj, where 
Ai = 60 MeV was obtained. 

We consider that the use of the w.f. at the origin |i?„(0)|gxp, extracted 
from di-electron widths, is the most promising because these w.f.s take im- 
plicitly into account the relativistic corrections as well as the influence of 
open channel(s), in this way drastically simplifying the theoretical analysis. 
A comparison of these w.f.s with those calculated in different models puts 
serious restrictions on the static potential used and also on many-channel 
models. 

We also show that the nonperturbative spin-spin interaction gives a con- 
tribution of about 9 MeV in the J/ip — 7]c{lS) mass difference. 



2 The spin-spin interaction 

The HF splitting between the n^Si and n^So levels will be considered here 
for two cases. The first one corresponds to the standard perturbative (P) 
spin-spin interaction with a 5-function: 

VLir) = s,.s, f^asifi) (l + ^p) 6ir) ^ s, ■ s.VuAr), (2) 

which in one-loop approximation gives the following HF splitting 0: 

KAnS) = ?^|fi„(0)r (l + ^ . (3) 

where p = -^Po — | — | In 2. This correction is small: < 0.5% in bottomonium 
{uf = 5) and < 3% in charmonium (jif = 4). 

It is very probable that S{r) may be considered as a limiting case and the 
"physical" spin-spin interaction is smeared with a still unknown "smearing" 
radius. For example, for the Gaussian smearing function, 

6{r)^^ fr^drexp{-(3^r^), (4) 



the splitting can be rewritten as 

^UnS) = l^UP)\Rnm' (l + -P) , (5) 

where by definition "the smearing factor" ^n{P) is 

UP) = -^tJ^ I \Rn{r)?eM-f^'r')r'dT. (6) 



The general expression Eq. (0) is evidently valid for other smearing pre- 
scriptions which may differ from Eq. (@)). The factor ^^ is calculated in 
Appendix ^ 

It is well known [H' that the w.f. at the origin is very sensitive to the 
form and parameters of the gluon-exchange interaction and also to the value 
of the quark mass used. Therefore we make the following remarks 



1. To miniinize the uncertainties in the w.f. at the origin, i?„(0), we shall 
use here the w.f.s extracted from the experimental data on the leptonic 
widths and denote them as |i?n(0)|exp- In this way the relativistic cor- 
rections to the w.f. and the influence of open channel(s) are implicitly 
taken into account. 

2. In Eqs. © and (0) the constituent mass Uq enters. This fact can be 
rigorously deduced from relativistic calculations [HI [T3] , 



UJr 



[nS) = {^p^ + mD^s, (7) 



where under the square-root the pole mass m^ = mg(pole) is present. 
This mass is known with good accuracy and we take here mf,{po\e) = 
4.8 ± 0.1 GeV and mc(pole) = 1.42 ± 0.03 GeV, which correspond to 
well-established current masses Thh{rhf,) = 4.3(1) GeV, rhcirhc) = 1-2(1) 
GeV [3 . The important feature of the constituent mass ujq{nS) is that 
it takes into account the relativistic correction and its values depends 
weakly on the quantum numbers and on the static interaction used. 
Indeed, they lie in a rather narrow range for all nS states, both in 
charmonium and bottomonium: 

uJc{lS) = 1.62 ± 0.03 GeV(r2 = 1), uJc{nS) = 1.71 ±0.03 GeV{n > 2), 
ujb{nS) = 5.05 ±0.15 GeV. (8) 

Note that just these mass values are mostly used in nonrelativistic 
calculations, thus in a hidden way taking into account relativistic cor- 
rections. 

3. The leptonic width of the n^Si state in heavy quarkonia is defined by 
the Van Royen-Weisskopf formula with QCD correction 7^ [TH] : 

ree(n3^l)|exp = -^ |i?n(0) |Lp7.. (9) 

n 

Here e^ = | (|j for a 6(c) quark; a = 1/137, M„ = M{n^Si), and the 
QCD factor is given by 

^q{nS) = 1 - — «.(2m,). (10) 



The renormalization scale /z in Eq. (jTUI) is supposed to be known, /i = 
2mq(pole), as in Refs. CD] and in rji, -^ 77 decay [T7j. In some cases 
/i = M„ is also taken. With an accuracy of < 1% in bottomonium and 
< 10% in charmonium both choices coincide and therefore one can take here 
/x = 2mg {2mb = 9.6 GeV, 2mc = 2.9 GeV). 

Since for nj = 5 the QCD constant ^jrq is well known from high en- 
ergy experiments [3], the factor ■jb is supposed to be determined with a 
good accuracy. For A^g--(3 — loop) = 210(10) MeV, which corresponds to 
aj^{Mz) = 0.1185, one finds 

lb = Ibn = 0.700(5), as{2mb) = 0.177(3). (11) 

In charmonium [uf = 4) the strong coupling ajYs{2mc = 2.9 GeV) is deter- 
mined with less accuracy and for A^^r-- = 0.260(10) MeV one obtains 

a,{2mc = 2.9 GeV) = 0.237(5), 7^ = 0.60(2). (12) 



Here, the theoretical error comes from the uncertainty in the value A^^^. 
Then the w.f. at the origin, extracted from the di-electron width in Eq. 






l^n(0)lLp=^^^^^S^^^^ (13) 



implicitly takes into account the relativistic corrections as well as the in- 
fluence of the open channels, which gives rise to smaller values for |i?n,(0)| 
and the HP splitting. An additional decrease of Ahf comes from a possible 
smearing of the S-wave spin-spin interaction (for the P-wave states this effect 
is very small). 

The extracted values of |i?„(0)|gxp in the bb and cc systems are presented 
in Tables HI and El 

With the use of these w.f.s at the origin and the values 7?, = 0.70 we also 
very precisely reproduce the ratios of the leptonic widths measured by the 
CLEO Collaboration (third Ref. in [12]): 

ree(T(2^))/r,,(T(15)) = 0.46(1), ree(T(3S))/ree(T(15)) = 0.32(1) 
ree(T(35))/ree(T(25)) = 0.69(2), (14) 

which confirms our correct choice of equal values of 7f,„ for these states. 

It is of interest to compare the extracted values of |-Rn(0)|exp with the the- 
oretically predicted values, which mostly depend on the strong coupling used 



Table 1: The w.f. \Rn{0)\l^p from Eq. ^ (in GeV^) and the leptonic widths 
ree(T(nS)) (in keV) for the T{nS) states (7^ = 0.70) '^) 





1 ee('^*-' jexp 


iRnmi,'^ 


IS 


1.314(29) 
1.336(28) 


7.094(16) 
7.213(15) 


2S 


0.576(24) 
0.616(19) 


3.49(15) 
3.73(12) 


3S 


0.413(10) 


2.67(7) 


4S 


0.25(3) 


1.69(20) 


5S 


0.31(7) 


2.21(49) 


6S 


0.13(3) 


0.95(22) 



") For the states IS* and 25" the upper entries are taken from PDG 0| and 
the lower ones from the CLEO data ^]. The numbers concerning the 35* 
state are also taken from the CLEO data [12]. The values for the states 45", 
55*, and 65* are taken from PDG j3]. 
^^ In the third column only experimental errors are given. 



Table 2: The w.f. |^n(0)|^^p 
(in keV ) in charmonium (7^ = 



;in GeV^) 
: 0.60)'^). 



and the leptonic widths Fg 



(n 



'S,) 





IS 


2S 


3S 


4S 


ree(exp) 


5.40(22) 
5.68(24) 


2.12(12) 
2.54(14) 


0.75(1) 
0.89(8) 


0.47(15)) 
0.71(10) 


^n(0)|Lp 


0.911(37) 
0.959(40) 


0.51(3) 
0.61(3) 


0.22(1) 
0.26(2) 


0.16(5) 
0.24(4) 



"^ The upper entries are taken from Ref. j2j and the lower ones are taken 
from Ref. [13]. 



in the gluon-exchange (GE) term. In particular, if the asymptotic freedom 
(AF) behavior of astatic (^) is neglected, then the theoretical values can be 
2 — 1.5 times larger than |_R„(0)|gxp, even for those T{nS){n = 1, 2, 3) states, 
which lie far below the BB threshold ^^ (see Table El in Appendix ^ . 

In our previous analysis of the spectra and fine structure splittings in 
heavy quarkonia j3|Tl|[Tn] we have used the static potential Vsir) in which 
the strong couphng in coordinate space asir) is defined as in background 
perturbation theory (see Eq. ()A.3|) ). In bottomonium (in the single-channel 
approximation) the potential VB{r) appears to give values of |-R„(0)|th very 



close to the numbers l-R^ 



2 

exp' 



^n 



For illustration in Table 01 the ratios 

|-Rn(0)|exp 



l^n(0)L\ 



(15) 



are given for all known nS levels in charmonium and bottomonium. 

As shown in TableEl in bottomonium for the potential Vsir) the influence 
of open channels appears to be important only for the 45* and 65* levels, while 
for other states the single-channel calculations are in good agreement with 



Table 3: The factor Sn defined by Eq. (jl5|) for the potential VB{r) from 
Eq. ()A.2|) for the T(nS) and ip{nS) wave functions at the origin"). 





IS 


2S 


3S 


4S 


5S 


6S 


hh 


1.08(4) 


1.02(4) 


1.02(4) 


0.72(9) 


1.03(22) 


0.47(10) 


cc 


1.01(4) 


0.82(5) 


0.41(2) 


0.32(10) 







"''' Given numbers correspond to Fee taken from PDG [S]. 

experiment. It is not so for many other potentials [H] and this means that 
any conclusions about the role of open channels cannot be separated from 
the qq interaction used in a given theoretical approach. 

In charmonium the effect of open channels is much stronger and reaches 
~ 60% for the 3^ and the 4^ states {Sn = 0.4) and about 20% for the 
i/j{2S) meson. This number is even smaller (12%) if the new CLEO data 
from Ref. [12] are used. 



3 The hyperfine splittings in bottomonium 

We consider two cases: 

A. No smearing effect, i.e. in Eq. (0) the smearing parameter ^bn = 1-0. (for 

all n). 

B. The smearing parameter ^b„, Eq. (P), is calculated with the value /3 = 

-\/l2 GeV, corresponding to the smearing radius r^s = /3"^ = 0.057 fm. 



Unfortunately, at present there is no precise prescription as to how to 
choose the renormalization scale in the HP splitting Eq. Q: in cxjjg{tl) the 
scale ji = mfe(pole) = 4.80 ± 0.01 GeV is often used. With Aj^{nf = 5) = 
210(10) MeV (just the same as in our calculations of 7;, Eq. (fTT]) ) we find 



as{bb,jl) = ajjs{4-8 GeV) = 0.21(1) 



(16) 



Table 4: A^p{nS) (in MeV) in bottomonium for ajjg{Mn) = 0.21; Ub = 5.10 



GeV and |i?n(0)|^^p from Tabled 





^b = 1.0 


^bn from Appendix |Bl 

(3 = Vl2 GeV 
rss = 0.057 fm 


IS 


51 (4) (4) 


28 (2) (3) 


2S 


25 (3) (2) 


12 (2) (1) 


3S 


22 (5) (2) 


10 (2) (1) 


4S 


12 (3) (1) 


5.1 (2) (1) 


5S 


16 (2) (1) 


6.4 (1) (1) 


6S 


7 (2) (1) 


2.7(1)(1) 



With this value of as{fl) and |i?n(0)|exp from Table HJ one obtains the HF 
splittings in bottomonium presented in Table |3J second column. 

The numbers in Table |3] contain an experimental error coming from 
TeelnS) [3] (first number), and a theoretical error (second number) < 10% 
( < 4% from the Ub value and < 5% comes from aj^sii^t))- For the T(nS) 
states (n = 1, 2, 3) the calculated HF splittings (^^ = 1.0) are very close to 
the splittings from Refs. jH]. 

In TablelHthe HF splittings, calculated with the smearing function Eq. (jH)), 
are also given (third column). The smearing radius, Tss = (3~^ = 0.057 fm 
(/? = v^l2 GeV), is taken to fit the experimental values of the J/-ip — r]c{lS) 



and ip{2S) —ric{2S) splittings (see the next section). However, even for such a 
small radius Ahf('^5') turns out to be 50% (n = 1 —4), 60% (n = 5, 6) smaller 
compared to the "nonsmearing" case. In particular, the T{lS) — rib{lS) split- 
ting turns out to be 28 MeV instead of 51(4) MeV for ^bn = 1-0. For higher 
excitations very small splittings, Ahf ~ 6 MeV and 3 MeV, for the 55* and 
65* states are obtained. 

Our value of r^g = 0.057 fm is very close to that from Ref. ^Ul where 
rss = 0.060 fm is taken, both in charmonium and in bottomonium. However, 
in spite of this coincidence our numbers are about two times smaller than in 
[TU], where T(IS') —r]b{lS) = 60 MeV is obtained. For the 25" state our value 
of splitting is 12 MeV, still smaller than the value 20 MeV found in Ref. pH] . 

We can conclude that observation of an ribijiS) meson could clarify the 
role of smearing in the spin-spin interaction between a heavy quark and 
antiquark. 

4 r]c(nS) masses 

The splitting Eq. (0) factually depends on the product as{fi) x ^„, therefore 
it is convenient to discuss an effective HF coupling: 

aupinS) = asifin) ^m, (17) 

which is the only unknown factor since |-R„(0)p can be extracted from the 
di-electron widths and the constituent masses uJc{nS) are known with < 5% 
accuracy from relativistic calculations. These masses may be specified for 
different nS states, Refs. [TH [T3] : 

io^ilS) = 1.62(3) GeV, uj^{2S) = 1.71(4) GeV, 
uj,{3S) ^ uj^{4S) = 1.73(4) GeV. (18) 

Here the theoretical errors come from a variation of the parameters of the 
static potential. 

As discussed in Ref. [7], the experimental splittings J/ip — ric{lS) and 
ip{2S) — ric{2S), can be described if a different auF Eq. (|T7jl for the IS" and 
2S states are taken, namely, aHF(lS) — 0.36 and aHF(2S) = 0.30. Such a 
choice implies two possibilities: 

A a^(/ii) = 0.36, as(/i2) = 0.30, ^^(/is) = asifJ'4.) < 0.30, ^^n = 1-0 
(for all n). i.e., the renormalization scale is supposed to grow for larger 

10 



excitations. In particular, for the value AJ^(2 — loop) = 270 MeV 
we have /xi = 1.25 GeV = rhc^rhc), while the scale fi2 = 1.60 GeV is 
essentially larger. For this choice of ohf the perturbative HF splittings 
are given in Table El 

B The normalization scale /i„ (and therefore the coupling constant as{fin)) 
is supposed to be equal for all nS'-states: as{fin) = 0.36. In this case 
a smearing of the spin-spin interaction is of principal importance to 
explain the experimental value for the ipi^S) — ric{2S) splitting. 

Besides, we have also calculated the contributions coming from the NP 
spin-spin interaction. In bottomonium their values are small, A^p(ra5') < 1 
MeV and can be neglected. In charmonium, as well as in light mesons, 
the situation is different. Due to the NP spin-spin interaction in the IP cc 
state a cancellation of the perturbative and NP terms takes place ^Hl- As a 
result, the mass difference Mcog(xcj) — M{hc) = (1 ± 1) MeV turns out to 
be close to zero or even positive, in accord with experiment ^Hl- Just the 
same NP contribution, Eq. (fT^. provides the correct value of the splitting 
^cog(aj) - M(6i(lP)) ^ 22 MeV in light mesons 120]. 

The values of A^p(nS'), (in Table El fourth column) are calculated in 
Appendix O They can be slightly different for different values of the gluonic 
correlation length Tg which defines the behavior of the vacuum correlation 
functions (v.c.f.). For Tg = 0.3 fm the v.c.f. has an exponential behavior 
over the whole range, < r < cx3 [21] and in this case 



.NP/ o\ " '-'2 



18cj2 



A^^(nS) = —^1.20(1 ± 0.07)J^(n^), (19) 



where the number 1.20(1±0.07) follows from lattice data |2I]. In Eq. (jl9|) the 
gluonic condensate G2 is taken to be equal to 0.043(3) GeV^, as in Ref. fIE\ . 
and the matrix element (m.e.) is defined as 

J^inS) = {rKi {r/Tg))ns. (20) 

Our calculations give 

T{nS) = 0.80 GeV-^(15), 0.40 GeV-\2S), 

0.27GeV-^(3S), 0.20 GeV-^(45). (21) 

The corresponding A2p(raS') are given in TableEl One can see that A^f(1>S') = 
9 ± 2 MeV needs to be taken into account in the J/ip — ric{lS) splitting while 

11 



Table 5: The splittings A|^p(n5') and AyF(?T-5') (in MeV) in charmonium"-' . 





APp(n5) 

(no smearing: ^c = 1-0) 

a,(/ii) = 0.36 

a,(/i„) = 0.30 

(n = 2,3,4) 


APp(n5) 

Tss = 0.29 GeV-i 

dsiHn) = 0.36 

(n = 1 - 4) 


G2 = 0.043 GeV^ 


15 


118(5) 


100(6) 
108(7)^=) 


9±2 


experiment 


117(2) 


117(2) 




25 


51(5) 
61(5)^=) 


46(3) 
55(4)^=) 


3.5 ±1.5 


experiment 

^(25) - 7],{2S) 


48(4) 


48(4) 




35 


21(2) 


16(2) 


2±1 


45 


15(4) 


12(4) 


1.5 ±0.5 



"■^ The w.f. |-Rn(0)|exp is taken from Tableland corresponds to Tef,{nS) from 

PDG 0. 

'') The NP splittings are calculated in Appendix O 

'^^ Here |^i(0)|2^p = 0.959 GeV^ and |^2(0)|^^p = 0.61GeV3 are extracted 

from the CLEO data ^. 
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for the higher states the values of A^p he within the experimental and theo- 
retical errors. Their values are Ag^ = 3(2) MeV, 2(1) MeV, and 1.5(5) MeV 
respectively for the 25, 3S, and 45* states. A smaller value of Tg{Tg < 0.2 
fm), which also cannot be excluded, practically does not change the num- 
bers. Adoption of this value would only slightly decrease those splittings. 
Thus one can conclude that in case A with different renormalization scales 
yUn(yUi ~ 1-25 GeV is small, /i2 = 1-6 GeV), for the J/ip — ric{lS) and ip 
(2S) —T]c{2S) splittings agreement with experimental values can easily be 
obtained. 

If the renormalization scales //„ are supposed to be equal for all nS 
states: as{fin — fric = 1-25 GeV ) =0.36, then to explain the relatively small 
ip{2S) — Tjci^S) splitting a smearing effect needs to be introduced. Then for 
the potential used the values in{,cc) = 0.85, 0.80, 0.78, 0.76 for the 15, 2S, 
3S, and 43 states respectively, are calculated in Appendix |Bl and for this 
case the values of A^p{cc,nS) are also given in Table El For the 3S and 4S 
levels the values we predict are about 21(15) MeV (no smearing) and 16(12) 
MeV (with smearing), i.e. the difference between the cases A and B is only 
~ 20%. Notice that in case B the NP contribution improves the agreement 
with experiment for the states J/ip — r]c{lS). 

Thus, in the ip{nS) — rjcijiS) splittings the smearing effect appears to be 
less prominent that in bottomonium. 

In Appendix ^ we also show that the relativistic correlations as well as 
the constituent masses stop to grow when the many-channel description is 
effectively used. This fact is very important for the study of higher excitations 
in charmonium. 

5 Conclusions 

In our paper we have shown that 

1. In bottomonium /S^-p{nS) appears to be very sensitive to smearing of 
the spin-spin interaction. Due to this effect the splitting decreases from 
51 MeV to 28 MeV for the 15 state and from 25 Mev to 12 MeV for 
the 25 state; very small values are obtained for higher states. 

2. In charmonium there are two possibilities to describe Ahf(15) and 
Ahf(25), which are known from experiment. The first one refers to a 
different choice of the renormalization scale: /ii =1.25 GeV and /i2 — 
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1.60 GeV for the IS" and 25" states, if the smearing effect is absent. The 
second possibihty imphes equal renormahzation scales /i„(?2 = 1 — 4) 
for all nS states. Then to explain the ip{2S) — rjc{2S) splitting the 
smearing of the spin-spin interaction needs to be taken into account. 
We also expect that for the 15 level a small contribution (~ 9 MeV) 
comes from the NP spin-spin interaction. 

3. The 'ip{'iS) — '/7c (35') splitting is predicted to be around 16(2) MeV, 
without and 12(4) MeV with smearing effect. 

To understand the true role of the smearing effect in the spin-spin interaction 
the observation of an rib{nS) is crucially important. 

A The wave function at the origin 

Two factors strongly affect the w.f. at the origin: 

(i) The AF behavior of the coupling in the GE term: Vo(r) = ar — ^ °''^'^ , 

(ii) The influence of open channels. 

For illustrations of these statements in Table IHl the w.f.s |i?„(0)|th are given 
for several cases. 

1. For the Cornell potential, Vc{r) = err — |^ in single-channel approxi- 
mation {ac = 0.39, a = (2.34)"^ GeV^,^^^ = 1.84 GeV). 

2. For the same Cornell potential while open channels are taken into 
account with the use of the nonrelativistic Cornell coupled- channel 
(CCC) model j22]- Then the w.f. in the coupled-channel system, 
tpni'^) = 'ipin^Si), is presented as the composition of different contribu- 
tions, e.g. 

^^{l^Si) = 0.983|15> -0.050|25> -0.009|35> -0.003|45>, 
^2(2^5i) = 0.103|15> +0.838|25> -0.085|35> 

-0.017|45 > -0.007|55 > +0.040|1D > -0.008|2L) >, 
^3(3351) = 0.02e-*°-°^"|15 > +0.19e-^°-=^°"|(25) > +0.67|35 > 

+0.07e^°-^^"|45 > +0.04e^°-^^"|lD > +0.04e^°-^^"|2D > . 

(A.l) 
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In single-channel approximation the w.f. |-R„5(0)p for the Cornell po- 
tential and the derivatives |-R^(0)p were calculated in jS]: l-Ri^lO)!^ = 
1.454 GeV^, \RgiO)\^ = 0.927 GeV^, \RgiO)\^ = 0.791 GeV^, \R'isW ■ 
0.725 GeV^; \R'Id{0)\1 = 0-030 GeV^; |i?2D(0)lc = 0-0655 GeV^ From 
Eq. (IXl|l it follows that in the CCC model the 2^51 - l^Di and the 
3^5*1 — 2^Di mixings appear to be small compared to the analysis in 
[231 5 where the admixture of the l^Di state is ~ 22%. 



3. For the potential 



3 r 
the coupling is defined as in Refs. P^ IT3]: 



„ smgr 



VB{r) = ar-^^^, (A.2) 



1 



with 

q' + Ml 



(3i \ntB 
Po ^B 



(A-3) 



tB{q) = In 



Alinj) 



Here Mb = 0.95 GeV, is the background mass, KBijif) is expressed through 
the AjY^{nf) and in 2-loop approximation Ab(^/ = 4) = 0.360(10) MeV P^ . 
The string tension cr = 0.18 GeV^ and the same values of Uf. = 1.65 GeV are 
taken here for simplicity for all nS states. 

In Table ini for comparison the values of |-Rn(0)|exp extracted from the 
di-electron widths are also given. Then one can see that in single-channel 
approximation very different |i?„(0)|^jj are obtained for K(r) and VB{r). For 
the potential Vb(?") the w.f. |i?i(0)p appears to be very close to the "exper- 
imental" number: |i?i(0)p = 0.91 for J/ip, while for the Cornell potential 
|i?n(0)|th is about two times larger than |-Rn(0)|exp (single-channel approxi- 
mation). Moreover, even in the CCC model |i?ri(0)|th is still too large, being 
~ 15%(30%) larger for the 2S'(1S') states compared to \Rn{0)\l^p. Corre- 
spondingly, for such multi-channel w.f.s the leptonic widths will be 15% and 
30% larger than the experimental value, respectively. 

At this point we would like to stress that one cannot fit the leptonic widths 
Tee{nS) taking a much smaller value for the factor 7c, Eq. (fTT|) (or larger as(/x) 
in 7c). Such a change in the renormalization scale would come in contradic- 
tion with the description of other "annihilation" decays like ric{nS) —>■ '-ff, 
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Table 6: Comparison of |-Rn(0)|exp (iii GeV'^) from Eq. (fT^ with the w.f. 
l-^n(0)|th for the Cornell potential (single-channel approximation and also in 
the coupled-channel model) and for Vsir) Eq. ()A.2|) . 







IS 


2S 


3S 


Cornell 
potential 


single-channel 
approximation 


1.454 


0.928 


0.790 


coupled-channel 
model [21] 


1.268 


0.703 


0.520 


Vsir), 
single-channel approx. 


0.900 


0.616 


0.534 


\Rn{0) 2^p from Table 2'^) 


0.91(4)exp(5)th 

0.96(4)exp(5)th 


0.51(3)exp(2)th 
0.61(3)exp(2)th 


0.22(3)exp(l)th 
0.26(3)exp(l)th 



"•^ See the footnote to Table |21 for the origin of these numbers. 
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where just the value as{^ = 2m^ = 0.24(1), as in our analysis is taken. 
Also we would like to notice that in the CCC model the 2^5"! — l^Di and 
3^ Si — 2^Di mixings are small, as seen from Eqs. ()A.1|) . and instead large 
admixtures to the 2S{3S) states come from the neighboring 15(25') states. 
For the Vb('") potential Eq. ()A.2|) we have obtained good agreement with 
the "experimental" |-Ri(0)|gxp (for the J/ip ~ 5% accuracy) and 20% dis- 
crepancy for -0(25). For higher states open channels give rise to a suppres- 
sion of |i?„(0)|^h (single-channel approximation) by ~ 20% for iIj{2S) and 
about 60% for the '?/'(4040) and ■?/'(4415) mesons. It is precisely the effect 
of the open channels that appears to be responsible for the screening of the 
GE interaction which was discussed in Ref. j21], where open channels have 
been considered effectively through flattening of the confining potential and 
switching off of the GE interaction. 

B The smearing factor ^^ 

For the smearing function Eq. (@)) the factor 

a = e(n5) = ^— ^, (B.l) 



is calculated here for the potential VBij) Eq. ()A.2|) . This factor appears to 
be weakly dependent on a variation of the mass iriq. In Table [3 its values are 
given for the parameter (3 = v^l2 GeV which corresponds to the smearing 
radius r^g = j3~^ = 0.057 GeV. For smaller (3 (larger Tss) the smearing effect 
is becoming even larger (the ^n are smaller). 

From Table [7| one can conclude that in bottomonium (due to smearing) 
the HF splittings decrease by a factor of two for all nS levels (ra = 1 — 6), 
while in charmonium the smearing effect is weaker and mainly important for 
the 35 and 45 states (see Tables El and Ej) . 

We would also like to make several remarks about relativistic corrections. 
Partly this correction is taken into account through the use of the constituent 
mass ujq - the average over the quark kinetic energy Eq. ((Zj). Actually, in the 
Hamiltonian, or the spinless Salpeter equation (SSE): 

(T + V)^P^^in) = Mo(nL)^^Jr) (B.2) 

we use the expansion of the square root near the point p^ + m'^ — uj^. Then 
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Table 7: The smearing factor ^„ Eq. (jHl) in bottomonium and charmonium 
for the potential Vsir) Eq. |0|) . a = 0.1826 GeV^, ^^(r) is defined by 
Eq. (EH). 





e„(66)») 


p = Vl2 GeV 


Uccr 
P = 2 GeV 


IS 


0.61 


0.85 


0.69 


2S 


0.55 


0.80 


0.62 


3S 


0.51 


0.78 


0.57 


4S 


0.49 


0.76 


053 


5S 


0.47 


0.74 


0.50 



'^) mfe = 5.1 GeV 
*) Uc = 1.70 GeV 



Table 8: The ratio r] = \R^si^)\'^/\^nsi^)\'^ ^^ charmonium for the hnear 
potential aor and for the flattening potential a{r) ■ r (ctq = 0.18 GeV^). The 
mass mc(pole) = 1.45 GeV. The parameters of cr(r) = cxofl ~ li^")] ^-^e taken 
fromRef. ^. 





IS 


2S 


3S 


4^ 


aor potential 


1.108 


1.196 


1.258 


1.310 


a{r) ■ r potential 


1.093 


1.153 


1.147 


1.045 



the kinetic term is 



Ti 



R 



p^ + nil 



2uj 



\ 



1 + 



p2 + rrfl 



^t 



oJi 



U}a + 



ni„ 



+ 



UJn 



P' 



U}„ 



= T, 



EA, 



(B.3) 



which is different from the standard nonrelativistic case. This expansion, or 
so-called Einbein approximation (EA), appears to be a good approximation 
to the SSE j2n]- While relativistic corrections are important, as in charmo- 
nium, it provides much better accuracy than the NR approximation, where 
the difference between the constituent mass ojq{nL) and niq =const. is ne- 
glected. Still, for the w.f. at the origin a small difference between the exact 
solution R^^{r) of the SSE and an approximate solution for the EA equation 
occurs 



{Tea + V)i)EA{r) = MEA{nL)i,EA{r 
As an example, in Table |H1 the ratio rj = \R^si^)\'^ 



pEA 



(B.4) 

(0)p are given for 

two confining potentials. The first one refers to the single-channel case with 
the string tension ctq = const. In this case the relativistic correction (the 
factor 1]) is growing for higher levels and reaches about 30% for the 45* state. 
However, if there are open channels, the linear potential is fiattening, 
and such a fiattening potential a{r) -r effectively takes into account the open 
channels. This feature provides an essential shift down for high excitations 
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[23 • Due to the flattening, already for the 35" state the factors rj as well 
as u!c{SS), stop to grow, i.e., an open channel affects the w.f. at the origin 
much more strongly than the relativistic correction. Since at present there is 
no relativistic multichannel theory and one cannot distinguish between both 
effects, we expect that rather small relativistic corrections (< 15%) will take 
place for all nS states in charmonium. Note that this relativistic correction 
does not affect the smearing factor C,n{cc). 

C Nonperturbative spin-spin interaction 

The NP contribution to the HF splitting Mcog(Xcj) — M{hc) has been dis- 
cussed in detail in Refs. [IHl 120] , where the following general expression has 
been obtained 

2 

Kf^P) = ?^4l-20(l ± 0.07)e^ T{nL). (C.l) 

Here the matrix element is given by 

J^inL) = {e{r - ro) r K^ {r/Tg))nL, (C.2) 

with Tg the gluonic correlation length, which is still not rigorously fixed in 
a lattice measurement (see the discussion in Ref. [^). The parameter tq in 
Eq. ()C.1|) characterizes the size of the region near r = where the vacuum 
correlation functions are not well defined. In lattice data for large Tg = 0.3 
fm {uf = 4) it is supposed that ro = and then the expressions ()C.1|) and 
(1021) go over into Eqs. (HHEOl)- If ro y^ 0, (e.g. for Tg = 0.2 fm when the 
value ro = Tg = 0.2 fm) an additional term needs to be added to A^p ()C.1|) . 
At small r, {r —>■ 0), the NP spin-spin potential was defined in Ref. |18j . 

VhT = ^1 ■ ^2FhT Vinr^O) = ^^G,{ro + Tg). (C.3) 

Then assuming that the v.c.f. has a kind of plateau near the origin, one ob- 
tains the following contribution from the region r < tq (Eq. (|T8|) in Ref. ^ 



TT 



2 



V^NP(r -. 0) = ^G,{ro + Tg)^rl - r^ 0(ro - r). (C.4) 
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Therefore, in the most general case the total NP contribution is 



vr^Go ( , _x 12. 



A^^(nL) = —^ 1.20(1 ±O.O7)e^«^(nL)0(r-ro) 



+2(ro + T,)(v/r2-r2)„^^(ro - r)| . (C.5) 



Note that for Tg = r^ = 0.2 fm the second term in Eq. ()C.5|) is important 
only for the S'-wave states, being small for the P-wave states, i.e. for the 
splitting Mcog(xcj) — M{hc) considered in Ref. [TH] . 
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